Scattering and transport statistics at the metal-insulator transition: A numerical 
study of the power-law banded random matrix model 



J. A. Mendez-Berimidez 

Instituto de Fisica, Benemerita Universidad Autonoma de Puebla, Apartado Postal J-48, Puebla 72570, Mexico 

Victor A. Gopar 

Depto de Fisica Teorica, Facultad de Ciencias, and Instituto de Biocomputacion y Fisica de Sistemas Complejos (BIFI), 
Universidad de Zaragoza, Pedro Cerbuna 12, E-50009 Zaragoza, Spain 



o 

> 

O 

O 

o 

m 

c 

CO 



T3 

C 

O 

Q : 

> 
00 

in 

(N 

o 
o 



Imre Varga 

Elmeleti Fizika Tanszek, Fizikai Intezet, Budapesti Muszaki es Gazdasdgtudomdnyi Egyetem, H-1521 Budapest, Hungary and 
Fachbereich Physik und Wissenschaftliches Zentrum fur Materialwissenschaften, 
Philipps Universitat Marburg, D-35032 Marburg, Germany 
(Dated: November 2, 2010) 

We study numerically scattering and transport statistical properties of the one-dimensional An- 
derson model at the metal-insulator transition described by the Power-law Banded Random Matrix 
(PBRM) model at criticality. Within a scattering approach to electronic transport, we concentrate 
on the case of a small number of single-channel attached leads. We observe a smooth crossover 
from localized to delocalized behavior in the average scattering matrix elements, the conductance 
probability distribution, the variance of the conductance, and the shot noise power by varying b (the 
effective bandwidth of the PBRM model) from small (6 <C 1) to large (b > 1) values. We contrast 
our results with analytic random matrix theory predictions which are expected to be recovered in the 
limit b —¥ oo. We also compare our results for the PBRM model with those for the three-dimensional 
(3D) Anderson model at criticality, finding that the PBRM model with b £ [0.2, 0.4] reproduces well 
the scattering and transport properties of the 3D Anderson model. 

PACS numbers: 03.65.Nk, 71.30.+h, 73.23.-b 



I. INTRODUCTION 

The study of systems at the Anderson metal-insulator 
transition (MIT) has been a subject of intensive research 
activity for several decades.— ~— In particular, much in- 
terest has been focused on the scattering properties of 
critical systems by analysing the probability distribution 
functions of the resonance widths T and Wigner delay 
times T w r-~— as well as the transmission or dimension- 
less conductance TJ£~— The distribution functions of T 
and r w have been shown to be related to the properties 
of the corresponding closed systems, i.e., the fractality 
of the eigenstates and the critical features of the MIT. 
On the other hand, at the MIT, the distribution of con- 
ductances w(T) has been found to be universal, i.e., size 
independent, but dependent on the adopted model, di- 
mensionality, symmetry, and even boundary conditions 
of the system. w(T) has been studied for systems in two 
and more dimensions with a large number of attached 
single-channel leads^ - — In fact, concerning the conduc- 
tance of one-dimensional (ID) systems and its statistical 
distribution, the regime of small number of leads has been 
left almost unexplored 

In the present work we study numerically several sta- 
tistical properties of the scattering matrix and the elec- 
tronic transport across disordered systems in one and 
three dimensions described by the Power-law Banded 
Random Matrix (PBRM) model at criticality and the 
three-dimensional (3D) Anderson model at the MIT, re- 



spectively. We stress that we concentrate on the case of a 
small number of attached leads each of them supporting 
one open channel. 

The organization of this paper is as follows. In the next 
subsection we define the PBRM model, the 3D Anderson 
model, as well as the scattering setup. We also define 
the scattering quantities under investigation and provide 
the corresponding analytical predictions from random 
scattering-matrix theory (RMT) for systems with time- 
reversal symmetry. These analytical results will be used 
as a reference along the paper. In Section II we analyse 
the average scattering matrix elements, the conductance 
probability distribution, the variance of the conductance, 
and the shot noise power for the PBRM model as a func- 
tion of its effective bandwidth b. In Section III we com- 
pare the results of the PBRM model at criticality with 
the scattering and transport properties of the 3D Ander- 
son model. Finally, Section IV is left for conclusions. 



A. The PBRM and the 3D Anderson models 

As we have mentioned above, in the present study we 
adopt two models, namely, the PBRM model at critical- 
ity and the 3D Anderson model at the MIT. 

The PBRM mode&SS describes ID samples of length 
L with random long-range hopping. This model is rep- 
resented by A x A (A = L) real symmetric matri- 
ces whose elements arc statistically independent random 
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variables drawn from a normal distribution with zero 
mean, (fly) = 0, and a variance decaying as a power 
law (|fly| 2 ) ~ (b/\i — j\) 2a , where b and a are parame- 
ters. There are two prescriptions for the variance of the 
PBRM model: the so-called non-periodic, 
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where the ID sample is in a line geometry; and the peri- 
odic, 
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where the sample is in a ring geometry. Field-theoretical 
consideration s 4 ' 30 ^ 1 and detailed numerical investiga- 
tions ^ 32 ' 33 have verified that the PBRM model undergoes 
a transition at a = 1 from localized states for a > 1 to 
delocalizcd states for a < 1. This transition shows all the 
key features of the Anderson MIT, including multifrac- 
tality of cigenfunctions and non-trivial spectral statistics 
at the critical point. Thus the PBRM model possesses 
a line of critical points b G (0, oo). We set a = 1 in our 
study, i.e., we work with the PBRM model at criticality. 

The 3D Anderson model with diagonal disorder is de- 
scribed by the tight-binding Hamiltonian (TBH) 



fl = £|n)VF n (n| + E l n >< m l 



(3) 



(n,m) 



where n = (n x ,n y ,n z ) labels all the N ~ L 3 sites of a 
cubic lattice with linear size L, while the second sum is 
taken over all nearest-neighbour pairs (n, m) on the lat- 
tice. The on-site potentials W n for 1 < n x ,n y ,n z < L 
are independent random variables. When W n are Gaus- 
sian distributed, with zero mean and variance W 2 /12, 
the MIT at energy E ~ occurs for W = W c ~ 21.3. 
See Refs. 011131. Then, for W < W c (W > W c ) 
the system is in the metallic (insulating) regime. We set 
W = W c in our study. 

We open the isolated samples, defined above by the 
PBRM model and the 3D Anderson model, by attach- 
ing 2M semi-infinite single channel leads. Each lead is 
described by the ID semi-infinite TBH 



fl lead = ^(|n)(n + l| + |n+l)(n|) 



(4) 



n=l 

Using standard methods 3 ^ one can write the scattering 
matrix (S-matrix) in the for m 11 ' 37 

S(E) = (^ T t ^ = 1 - 2* sin(fc) W T (E - Hcs^W , 

(5) 

where 1 is the 2M x 2M unit matrix, k = arccos(£?/2) 
is the wave vector supported in the leads, and Hog is an 
effective non-hcrmitian Hamiltonian given by 



H cB = H - e lk WW T 



(6) 



Here, W is an TV x 2M matrix that specifies the posi- 
tions of the attached leads to the sample. Its elements 
are equal to zero or e, where e is the coupling strength. 
Moreover, assuming that the wave vector k do not change 
significantly in the centre of the band, we set E = and 
neglect the energy dependence of H c ff and S. 



B. R.MT predictions for the Circular Orthogonal 
Ensemble 



Notice from Eqs. (Q][5]) that in the limit 6 ^ co the 
PBRM model reproduces the Gaussian Orthogonal En- 
semble. Therefore, in that limit we expect the statistics 
of the scattering matrix, Eq. ([5]), to be determined by 
the Circular Orthogonal Ensemble (COE) which is the 
appropriate ensemble for systems with time reversal sym- 
metry. Thus, below, we provide the statistical results for 
the average S'-matrix and the transport quantities to be 
analysed in the following sections, assuming the orthog- 
onal symmetry. In all cases, we also assume the absence 
of direct processes, i.e., (S) = 0. 

We start with the average of the S'-matrix elements. It 
is known that 



(\Saa'\ 2 ), 
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where (•) means ensemble average over the COE. 

Within a scattering approach to the electronic trans- 
port, once the scattering matrix is known one can com- 
pute the dimensionlcss conductance T = Tr(ttT) and its 
distribution w(T). For M = 1, i.e., considering two 
single-channels leads attached to the sample, w(T) is 
given by 
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< T < 1 

1 < T < 2 
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For arbitrary M, the predictions for the average value of 
T and its variance are 
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respectively. For the derivation of the expressions above 
see for example Ref. (38j . Another transport quantity of 
interest is the shot noise power P = (Tr(ttf —ttHtf)), 
which as a function of ftl reads^r— 
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FIG. 1: (Color online) Scattering setup. 2M leads, shown as 
blue (gray) lines, are attached to 2M sites (black dots) of a 
ID sample described by the (a) non-periodic and (b) periodic 
versions of the PBRM model. The case M = 2 with L = 11 
is shown as example. 



In the following sections we focus on (l^wl 2 ), (T), 
var(T), and P for the PBRM model and the 3D Anderson 
model, both at the MIT. In all cases we set the coupling 
strength e such that (S) ~ in order to compare our 
results, in the proper limits, with the RMT predictions 
introduced above. 



II. PBRM MODEL 

We attach the 2M leads to the first 2M sites of the 
ID sample described by the PBRM model. That is, in 
the non-periodic version of the PBRM model, Eq. (pj, 
we attach the leads at the boundary of the system. Sec 
Fig. nja). While in the periodic version, Eq. @, we 
attach them to the bulk. See Fig.fTJb). In the latter case, 
finite size effects are considerably reduced. However, the 
quantities we analyse below are L-independcnt once L is 
much larger than the number of attached leads for both 
versions of the PBRM model. 

We point out that our setup is significantly different 
from the one used in Refs. (28ll29| . where the conductance 
has also been studied using the PBRM model. There, 
for example, the leads in the M = 1 case are attached to 
sites which are separated a distance of L/2 and L in the 
periodic and non-periodic versions of the PBRM model, 
respectively. In such situation the scattering quantities 
are strongly L-dependent. 



A. Average scattering matrix elements 

First we consider the case M = 1, where the 5-matrix 
is a 2 x 2 matrix. In Fig. dja) we plot the ensemble 
average of the elements |Sii| 2 and |Si2| 2 as a function 
of the bandwidth parameter b, (|Sii| 2 ) b and (|<Si2| 2 ) b , 
for the periodic (black symbols) and the non-periodic 
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FIG. 2: (Color online) Black [red (gray)] symbols: Average 
S'-matrix elements (|Sii| 2 ^ and ^|Si2| 2 ) for the periodic [non- 
periodic] PBRM model at criticality as a function of (a) b 
and (b) D2 for M = 1. The blue (gray) dot-dashed lines 
correspond to 2/3 and 1/3; the RMT prediction for (IS11I 2 ) 
and (|5i2| 2 ), respectively, given by Eq. (0. The black [red 
(gray)] dashed lines in (a) are Eqs. (|13[) and (|14[) with S w 
2.5 [S « 2.2]. The black [red (gray)] dashed lines in (b) are 
Eqs. (Jl8j) and (JTSJ with S « 2.5 [8 « 2.2]. Error bars are not 
shown since they are much smaller than symbol size. 



(red symbols) PBRM model. We concentrate on these 
two matrix elements since the other two elements give 
no additional information: (\S22\ 2 ) b = (I'S'iil 2 ^ an d 
(\S2i\ 2 ) b = (|<5i2| 2 ) b - Notice a strong 6-dependence 
of the average S'-matrix elements driving them from a 
localized-like regime [(l-Siil 2 ) ps 1 and <( 1 5"i2 1 2 ) ~ 0; i.e. 
the average conductance is close to zero] for b <C 1, to 
a dclocalizcd-like or ballistic- like regime [(|£ii| 2 ) ~ 2/3 
and (IS12I 2 ) ~ 1/3; i.e., RMT results are already recov- 
ered] for b > 4. 

Moreover, we have found that (|Sii| 2 ) and (IS12I 2 ), as 
a function of the bandwidth 6, are well described by 



(|5 U | 2 > 6 = l-(|Si 2 | 2 ) b , 



<|Sd 2 ) fc = i 



1 



1 + (6b)- 



(13) 
(14) 



where 8 is a fitting parameter. Eq. (fTB")) is a consequence 
of the unitarity of the scattering matrix, SS^ = 1, while 
the factor 1 /3 in Eq. JT3) comes from Eq. with M = 1. 
In Fig. [2£a) we plot Eqs. (JT3) and COP (dashed lines) 
and compare them with the corresponding results from 
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the PBRM model (symbols) in the periodic and non- 
periodic setups. In the same figure Eq. (O is also plotted 
(dot-dashed lines). 

On the other hand, it is well known that in systems 
at the disorder driven MIT both the energy spectra and 
the eigenstates exhibit multifractal characteristics^ The 
PBRM model is characterized by the effective bandwidth 
b that drives the system from strong (b — > 0) to weak 
(b — > oo) multifractality. Multifractality can be quan- 
tified by the generalized dimensions D q which describe 
the fluctuations of the eigenfunctions. The multifractal 
dimensions D q of the cr-th eigenfunction \E ,cr (given as a 
linear combination of the basis states in a system with lin- 
ear size L, <3> CT = J2t=i M are defined through the so- 

,2g 



called inverse participation numbers, I q 
by the scaling"^ 



Er=i ic? 



(15) 



where (I q ) typ = exp(lnl g ) is the typical value of I q . 
However, among all dimensions, the correlation dimen- 
sion _D 2 plays a prominent role4^ As b transits from zero 
to infinity, D2 takes values from zero to one. 

Using numerical diagonalization of matrices with sizes 
L = 128, 256, 512, and 1024, we extracted D 2 by the 
use of Eq. ([To]). For each system size L, we used a num- 
ber of disorder realizations in order to have at least 10 5 
data for its analysis. For each disorder realization we 
used 25% of the states at the centre of the spectral band 
(where the density of states is almost constant) in or- 
der to avoid boundary effects. We found good agreement 
between the numerically obtained D2 and the analytical 
estimations^ 



2b , 

1 - (lib)- 1 



6<1 
6> 1 



(16) 



Moreover, in Ref. [44j the following phenomenological an- 
alytical expression for D2 as a function of b was proposed: 




ln(i7T ) 



FIG. 3: (Color online) (a) Black [red (gray)] curves: Probabil- 
ity distribution p(lnT) for the periodic [non-periodic] PBRM 
model at criticality for several values of b < 1 (b = 0.01, 
0.02, 0.04, 0.1, 0.2, and 0.4 from left to right) in the case 
M = 1. (b) Black [red (gray)] symbols: (lnT) as a function 
of b for the periodic [non-periodic] PBRM model. The blue 
(gray) dashed line is the best fit of the data to the logarith- 
mic function A + lnb 2 with A « -1.44. (c) p(lnT) for b < 1 
scaled to T typ = exp(lnT) ~ b 2 . The blue (gray) dashed line 
proportional to (T/Tty P ) is plotted to guide the eye. 



D 2 (b) 



1 



1 + (76)- 



(17) 



where 7 is a fitting parameter. See also Ref. [Tj]. Eq. p7|) 
describes well our numerical results for D 2 with 7 = 2.94 
and 7 = 2.88 for the periodic and non-periodic versions 
of the PBRM model, respectively. 

Then, by substituting Eq. dTTJ) into Eqs. (US]) and (fT4|) 
we obtain the following expressions for the averages of 
|5ii| 2 and 1 612 1 2 as a function of 



<|Snr% 2 = i-(|s 12 | 2 ) D2 , 
<l^l 2 k = \ — 
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(18) 
(19) 



is interesting to note that 



D 2 



D 2 ^0 
D 2 -> 1 



(20) 

which might be relevant for systems at the MIT where 
D2 can be tuned. 

Finally we want to remark that concerning ^|<S aa ' | 2 ) for 
the PBRM model, the RMT limit, expected for b ->• 00, 
is already recovered for b > 4. See Fig. [U^a). 



B. Conductance probability distribution 



In Fig.[2jb) we compare Eqs. (fT8]l and (p~9]> (dashed lines) 
with the ensemble-average numerical results at different 
values of Z?2- The agreement is excellent. In addition, it 



Now we turn to the conductance statistics. For b <C 1 
the conductance distribution w(T) is highly concentrated 
close to T = 0. So it is more convenient to analyse the 




FIG. 4: (Color online) Black [red (gray)] curves: Conductance 
probability distribution w(T) for the periodic [non-periodic] 
PBRM model at criticality for some large values of b in the 
case M = 1. Blue (gray) dashed lines are w(T) C oe; the RMT 
prediction for w(T) given by Eq. ([5|. 

distribution of the transmission logarithm, p(lnT), in- 
stead. Then, in Fig.[3ta) we show p(lnT) for several val- 
ues of b < 1 in the case M = 1, for the periodic and non- 
periodic versions of the PBRM model. Notice that the 
distribution functions p(lnT) do not change their shape 
by increasing b, mainly for b <C 1, thus being scale in- 
variant. In fact, (InT) for b < 1 clearly displays a linear 
behavior when plotted as a function of In b as shown in 
Fig. [3{b). Then, all distributions functions p(lnT) fall 
one on top of the other when shifting them along the 
x-axis by the typical value of T, 

T typ =exp(lnT) , (21) 

as shown in Fig. GUJc). 

Note that for T < T typ , p(lnT) is proportional to 



L 


b = 0.1 


b = 1 


b = 10 


50 


-5.887 ±0.002 


-2.249 ± 0.002 


-2.001 ± 0.002 


100 


-5.886 ±0.002 


-2.247 ± 0.002 


-1.999 ±0.002 


200 


-5.875 ±0.002 


-2.245 ± 0.002 


-1.999 ±0.002 


400 


-5.875 ±0.002 


-2.244 ± 0.002 


-2.002 ± 0.002 


800 


-5.869 ±0.002 


-2.228 ± 0.002 


-1.989 ±0.002 



TABLE I: (In T) for the periodic PBRM model at criticality 
for some values of b and L in the case M = 1. 
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FIG. 5: (Color online) (a) Black [red (gray)] curves: Probabil- 
ity distribution p(lnT) for the periodic [non-periodic] PBRM 
model at criticality for several values of b < 1 (b = 0.01, 
0.02, 0.04, 0.1, and 0.2 from left to right) in the case M = 2. 
(b) Black [red (gray)] symbols: (InT) as a function of b for 
the periodic [non-periodic] PBRM model. The blue (gray) 
dashed line is the best fit of the data to the logarithmic func- 
tion A + lnb 2 with A w -0.057. (c) p(lnT) for b < 1 scaled 
to T typ = exp(lnT) ~ b 2 . The blue (gray) dashed line pro- 
portional to (T/T ty p) 2 is plotted to guide the eye. 



(T/T^p) 1 / 2 , see Figs.EJa) andEJc). Moreover, we found 
that such behavior extends from small to large values of 6, 
b 1. Also, notice that the behavior p(lnT) cx T 1 / 2 for 
small T coincides with the RMT prediction w(T) 
T -1 / 2 , sec Eq. (|HJ), since the change of variable T —> InT 
leads to p(lnT) = Tw COE (T). 

In Fig. |1 we show w(T) for large b (b > 0.4). In the 
limit b — > oo, w(T) is expected to approach the RMT 
prediction of Eq. ([8]). However, once b > 4, w(T) is 
already well described by w(T) COE . 

We want to stress that our results do not depend on 
the system size L, as can be seen in Table U where we 
report (In T) for the periodic PBRM model at criticality 
for some values of b and L in the case M — 1. In fact, 
(InT) does not depend on L out of criticality (a > 1 
and a < 1) either. This makes an important difference 
with respect to the setup where the leads are attached 
to opposite sides of the sample where (InT) increases 
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FIG. 6: (Color online) Black [red (gray)] curves: Conductance 
probability distribution w(T) for the periodic [non-periodic] 
PBRM model at criticality for some large values of b in the 
case M = 2. Blue (gray) dashed lines are w(T) G oe', the RMT 
prediction for w(T) given by Eq. ©. 



(decreases) as a function of L for a < 1 (a > 1 and 
a = 1). 

Now, in Figs. [5] and |H] we explore w(T) in the case 
M = 2. As well as in the case M = 1, studied above, 
here: (i) for small b, p(lnT) is scale invariant with T typ 
as scaling factor, sec Fig. [5j (ii) for b > 4, w(T) is well 
described by Eq. (0) , the corresponding RMT prediction, 
see Fig. [5] However, although p(\nT) for small b and 
w(T) for large 6 are practically the same for the periodic 
and non-periodic versions of the PBRM model, they show 
differences for intermediate values, 0.1 < b < 4, as can 
be seen in Fig. [SJa) and the upper panels of Fig. [5] 

We point out that for small T, T < T typ , p(lnT) is 
proportional to (T/T typ ) 2 , as shown in Fig. [5jc). This 
behavior is universal for the PBRM model with M = 2; 
i.e. it is 6-independent and valid for the periodic and non- 
periodic versions of the model. Again, as for the M = 1 
case, here for M = 2, the dependence p(lnT) oc T 2 in 
the limit b — > oo is consistent with the RMT prediction 
w(T) COE oc T for T < 1; sec Eq. ©. 

From the results above and since^ 



w(T) COE oc 7^ /2 /2-i 



(22) 



in the region < T < 1, we conclude that for the PBRM 
model 
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FIG. 7: (Color online) (a) Average conductance (T) as a func- 
tion of M for the non-periodic PBRM model at criticality for 
several values of b. (T) for the 3D Anderson model at critical- 
ity (3D AM) is also shown, (b) (T) / {T) COE as a function of 
b8 for the periodic PBRM model at criticality for M £ [1,5]. 
Insert: S versus M. S is obtained from the fitting of Eq. (|24[) 
to the (T) vs. b data. Thick full lines correspond to (T) = 0. 
Blue (gray) dashed lines are (a) the RMT prediction for (T) 
given by Eq. (fT0)l : and (b) one. 



for T <C T typ . We argue that the behavior dictated by 



pQnT) ocT A/2 / 2 



(23) 



Eq. (|23|) can be interpreted in the following way: Since 
small values of conductance T< 1 mean very strong re- 
flection, the scattering process is almost direct, i.e. an 
incident electron is scattered out the system mainly by 
the interaction with the sites at which the leads are at- 
tached. Then the electron does not explore the complete 
sample, not even part of it, and as a consequence it does 
not realize that the sample is at criticality. That is, for 
T <C 1, the incident electron docs not distinguish between 
a critical system and a random system represented by a 
full random matrix if the number of attached leads 2M is 
not much larger than 6; i.e., when the leads are attached 
to sites interconnected by Hamiltonian matrix elements 
located within the bandwidth of the matrix where all el- 
ements have almost the same variance: (|#ij| 2 ) ~ 1/2. 
So, we can use Eq. (|2"2")l and as a consequence Eq. (f2"3")) 
too. In fact, we have verified that Eq. ([23)) is valid up to 
M = 2 for b < 1, while for b = 10 it describes perfectly 
the left tail of p(lnT) for all the values of M used in this 
work (up to M = 5). Our result, given by Eq. (|2"3"]l . is 
the generalization of the result shown in Ref. [28[ , where 
the behavior p{hxT) oc T 1 / 2 was reported for the PBRM 
model in the localized (a > 1) and delocalized regimes 
(a < 1) for M = 1. 
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FIG. 8: (Color online) (a) Variance of T as a function of M for 
the non-periodic PBRM model at criticality for several values 
of b. var(T) for the 3D Anderson model at criticality (3D AM) 
is also shown, (b) var(T)/var(T) C oE as a function of b8 for 
the periodic PBRM model at criticality for M 6 [1, 5]. Insert: 
S versus M. 8 is obtained from the fitting of Eq. Q24p to the 
var(T) vs. b data. Thick full lines correspond to var(T) = 
0. Blue (gray) dashed lines are (a) the RMT prediction for 
var(T) given by Eq. and (b) one. 



FIG. 9: (Color online) (a) Shot noise power P as a function 
of M for the non-periodic PBRM model at criticality for sev- 
eral values of b. P for the 3D Anderson model at criticality 
(3D AM) is also shown, (b) P/Pcoe as a function of bS for the 
periodic PBRM model at criticality for M € [1,5]. Insert: S 
versus M. 5 is obtained from the fitting of Eq. l)24p to the 
P vs. b data. Thick full lines correspond to P = 0. Blue 
(gray) dashed lines are (a) the RMT prediction for P given 
by Eq. (TTJJ ; and (b) one. 



C. Mean and variance of the conductance and the 
shot noise power 

We now increase further the number of attached 
leads4£ In Figs. [7(a), E[a), and M, a ) we plot the aver- 
age conductance (T), the variance of the conductance 
var(T), and the shot noise power P for the non-periodic 
version of the PBRM model for several values of b with 
ME [1, 5]4£ It is clear from these three plots that chang- 
ing b from small (b <C 1) to large (b > 1) values produces 
a transition from localized- to delocalized-like behavior 
in the scattering properties of the PBRM model. That 
is, (i) for b < 0.1, (T) w 0, var(T) w 0, and P w 0; and 
(ii) for b > 10, (T). var(T), and P are well described by 
the corresponding RMT predictions given by Eqs. (|T0l) . 
(fTTj) . and (fF2]) , respectively. Similar plots are obtained 
(not shown here) for the periodic PBRM model. 

Moreover, we have observed that (T), var(T), and P 
behave (for all M) as {\Si 2 \ 2 ) b does. See Eq. Thus, 
we can write 





10 


p 




a 


10 








10 






10 



M= 1 



M = 2 



-,1/2 



-15 



-10 

InT 



FIG. 10: Probability distribution p(lnT) for the 3D Anderson 
model at criticality in the cases M = 1 and M = 2. Dashed 
lines with slopes T 1//2 and T 2 are plotted to guide the eye. 



average values, now for the periodic PBRM model as a 
function of b8 for M g [1,5]. Also, similar plots are 
obtained (not shown here) for the non-periodic PBRM 
model. 



X(b) = X COE 



1 



1 + (56)- 



(24) 



where X represents (T), var(T), or P and 6 is a fitting 
parameter. Then, in Figs. [TJb), [5(b), and [5(b) we plot 
(T), var(T), and P normalized to their respective COE 



III. 3D ANDERSON MODEL 

Since the most prominent realization of systems that 
undergo a MIT is the 3D Anderson model, it is of rele- 
vance to analyse its scattering and transport properties 
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FIG. 11: (Color online) Conductance probability distribution 
u>(T) for the 3D Anderson model at criticality in the cases 
(a) M = 1 and (b) M = 2. The red (gray) dashed line 
corresponds to w(T) for the non-periodic PBRM model with 
b = 0.2. 



taking as a reference the results shown in the previous 
section for the PBRM model. 

We attach the 2M leads to 2M sites at one of the edges 
of the cubic lattice described by the 3D Anderson model. 
In this way we make a line contact as we did in the case 
of the PBRM model. So, we can compare the scattering 
properties of both models. We use Gaussian distributed 
on-site potentials and system sizes from L = 6 to 10 (we 
have verified that our results do not change by increasing 
L further) with 10 6 to 10 4 disorder realizations. 

We start our analysis by looking at the average scat- 
tering matrix elements for M = 1. For the 3D Ander- 
son model at criticality we found (|5n| 2 ) « 0.926 and 
(IS12I 2 ) « 0.074, see TableOH which are close to those of 
the non-periodic PBRM model with b = 0.24. 

In Fig. [TU] we show conductance probability distribu- 
tions /5(lnT) for the 3D Anderson model at criticality in 
the cases M = 1 and M = 2. We found (InT) « -3.47 
and —3.06 for M = 1 and M = 2, respectively See Ta- 
ble HTH These values of (InT) for the 3D Anderson model 
are close to those of the PBRM model with b = 0.36 and 
b = 0.22, respectively. Notice also that, as well as for 
the PBRM model, for the 3D Anderson model p(lnT) in 
the region of T <C T typ is proportional to (T /T^) 1 / 2 and 
(T/T typ ) 2 for M = 1 and M = 2, respectively (dashed 
lines in Fig. ITUl). 

Additionally, in Fig.[TT]wc plot w(T) for the 3D Ander- 
son model in the cases M = 1 and M = 2. In the same 
figure we have also plotted (dashed lines) the conduc- 
tance distributions from the PBRM model with b = 0.2. 
We can see that the conductance distributions of both 
models are similar at this bandwidth value. 

We also compute (T), var(T), and P for the 3D Ander- 
son model and plot them in Figs.[7]Ja),[8][a), and[9ta), re- 
spectively (red curves labeled as 3DAM) . We observe that 
for the 3D Anderson model (T), var(T), and P behave as 
the corresponding quantities for the PBRM model with 
b close to 0.2. 

We want to stress that our results for the 3D Ander- 
son model at criticality do not seem to depend on the 
on-site potential distribution. Here we used Gaussian 



L 


/ 1 i2\ 


1 \ a i2\ 


6 


0.9256 ±0.0001 


0.0743 ± 0.0001 


8 


0.9254 ± 0.0003 


0.0745 ± 0.0003 


10 


0.926 ± 0.001 


0.073 ± 0.001 



TABLE II: Average S-matrix elements (ISii! 2 ) and (|Si 2 | 2 ) 
for the 3D Anderson model at criticality for some values of L 
in the case M = 1. 



L 


(In T) (M = 1) 


(InT) (M = 2) 


6 


-3.471 ±0.001 


-3.060 ± 0.001 


8 


-3.483 ± 0.005 


-3.058 ± 0.003 


10 


-3.450 ±0.016 


-3.078 ±0.012 



TABLE III: (InT) for the 3D Anderson model at criticality 
for some values of L in the cases M = 1 and M = 2. 

distributed potentials. However, we observe practically 
the same results by the use of box distributed on-site 
potentials (not shown here). 

IV. CONCLUSIONS 

We study the scattering and transport properties of 
the Power-law Banded Random Matrix (PBRM) model 
and of the three-dimensional (3D) Anderson model, both 
at criticality. 

We observed a smooth crossover from localized- to 
delocalized-like (ballistic-like) behavior in the scattering 
properties of the PBRM model by varying b from small 
(6 <C 1) to large (b > 1) values. For this crossover region 
we proposed heuristic analytical expressions for (\S aa ' | 2 )& 
and (|5 aa ' | 2 ) d 2 ■ F° r small b, we have shown that p(lnT) 
is scale invariant with the typical value of T, T typ , as 
scaling factor. We realized that RMT results, expected 
in the limit b — > 00, are already recovered for relatively 
small values of the bandwidth: b > 10. However, this fact 
is closely related to the small number of attached leads 
we used in this work: the larger the number of leads the 
larger the value of b needed to approach RMT behavior. 

Our conclusions are valid for leads attached to the bulk 
of the system as well as for leads attached at the bound- 
ary. In this work we assumed that time-reversal symme- 
try (f3 = 1) is present in our disordered systems. More- 
over, the case of broken time-reversal symmetry (/3 = 2) 
has been preliminarily studied in Ref. [45j where similar 
conclusions to this work were made. 

We have also shown that the scattering properties 
of the 3D Anderson model are similar to those for the 
PBRM model with b € [0.2,0.4]. This result is in agree- 
ment with previous studies^ 2 - were it was shown that sev- 
eral critical quantities related to the spectrum and eigen- 
states of the PBRM model with b = 0.3 are practically 
the same as for the 3D Anderson model. This makes the 
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PBRM model an excellent candidate to explore the prop- 
erties of the 3D Anderson model at criticality at a low 
computational cost. 

We want to stress that for both models at criticality, 
the PBRM model and the 3D Anderson model, we found 
the universal behavior pQnT) cx T M ' !/2 for T < 1, that 
we tested here for M = 1 and M = 2. Moreover by the 
use of the PBRM model with j3 = 2 (not shown here) 
we have already verified that the more general expres- 
sion p(lnT) cx T^ 2 / 2 , derived frorn^ W (T) cx T^ 2 / 2 " 1 
with < T < 1, holds. 

We emphasize that, even though we used a scattering 
setup where the leads are attached perfectly to sites at 
one side of the sample, our conclusions are not restricted 
to this topology4£ 

Finally, we recall that we concentrate here on the case 
of a small number of single-channel attached leads (up 
to ten). The study of the scattering and transport prop- 



erties of the PBRM and the 3D Anderson models in the 
regime of M ^> 1 will be the subject of a forthcoming 
contribution^ 
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All quantities reported in Figs. [2] to [6] were obtained for 
L = 50 using 10 6 ensemble realizations. In Subsection IIC 
we used L = 200 and 10 5 ensemble realizations. We have 
verified that in both cases our results are invariant by fur- 
ther increasing L. See also Table Q] 

We recall that for M = 5, ten single-channel leads are 
attached to the sample. 

In Ref. [2^] some properties of the transmission for the 
PBRM model were already studied for M 1. 



